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On the evolutionary stage of the unclassified B[e] star CD-42°11721* 
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ABSTRACT 

The star CD-42° 1 1721 is a curious B[e] star sometimes pointed as an evolved B[e] super- 
giant and sometimes as a young HAeBe star, due to very uncertain or even unknown stellar 
parameters, especially the distance. In this paper, we present new data gained from high- 
resolution optical spectroscopy and a detailed description of IR data of this star. We present a 
qualitative study of the numerous emission lines in our optical spectra and the classification 
of their line profiles, which indicate a non-spherically symmetric circumstellar environment. 
The first real detection of numerous [Fell] emission lines and of many other permitted and 
forbidden emission lines is reported. From our optical data, we derive an effective tempera- 
ture of T eS = 14 000 ± 1 000 K, a radius of i?* = 17.3 ± 0.6 R , as well as a luminosity of 
= (1.0±0.3) x 10 4 L Q . We advocate that CD-42° 11721 might be a post-main sequence 
object, even though a pre-main sequence nature cannot be ruled out due to the uncertain dis- 
tance. We further found that the SED in the optical and infrared can best be fitted with an 
outflowing disk-forming wind scenario rather than with a spherical symmetric envelope or 
with a flared disc, supporting our tentative classification as a B[e] supergiant. 

Key words: line: identification- stars: emission-line, Be - stars: individual: CD-42°11721. 



1 INTRODUCTION 

CD-42° 11721 (V921 Sco, Hen 3-1300, IRAS 16555-4237) is a 
very interesting galactic southern object showing the B[e] phe- 
nomenon. B[e] stars are stars of spectral type B with an optical 
spectrum showing strong Balmer line emission in combination with 
a large amount of permitted Fe II lines and forbidden O I and Fe II 
lines. In addition, these stars exhibit a strong near- and mid-IR ex- 
cess due to circumstellar dust. 

This star was observed by Merrill & Burwell (1949), who 
noted a peculiar natu re in its spectrum. Due to the presence of a 
nebulosity around it (Glass & Allen 1975) a nd possible spectral and 
photometric variations ( iHerbstl 1 1 9751 ; ICarlson & Henize|[l979h . it 
was tentatively classified as an HAeBe star (de Winter & The 1990; 
The, de Winter & P erez 1994). B ased on a possible high tempera- 
ture and luminositv. lVoorsI d 19991) has suggested that CD-42° 1 1721 
might be a young star but not a pre-main sequence object, 
while spectral similarities among CD-42 1 1721 an d CPD-52°9243 
(a supergiant candidate, Carl son & He niz el 1 19791) and HD 87643 
(an unclassified B[e] star, Huts emekers & Van Dronj Il990l) in- 
creased the confusion about the evolutionary state of this star. 
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An ev olved nature was reinforced by Hutseme kers & Van Droml 
(1990), who claimed the existence of a nebular overabundance 
of N/O. Other indications for a post-main sequence evolution- 
ary ph ase came from the proposed high lu minosity of this ob - 
ject bv lMcGregor, Hvland & Hillierl ( fl988l) and lshore et al.l l fl990h . 
This zoo of proposed evolutionary phases for CD-42° 11721 has 
lead to the i nclusion of this sta r into the list of uncl a ssified 
stars bvlLamers et al.l dl998l). Recently lHabart et al.1 d2003l) 
Hamaguchi, Yamauchi & Koyama (2005) suggested that CD- 



B[e 
and 

42° 11721 might be member of a cluster of several low and inter- 
mediate mass sources. 

The confusion concerning the evolutionary stage of CD- 
42 11721 is strictly linked to the absence of reliable physical pa- 
rameters. Since the optical spectrum does not display any photo- 
spheric lines, the spectral type derived using different sets of obser- 
vations and methods, varies between B0 and Aep, and consequently 
the T f ff ranges from 31600K to 12 3 00 K dHillenbrand et al.l 
Il992l ; ICidale. Zorec & Tringanielldl20"0lh . The range in distance 
is even worse: we found values from 1 36 pc up to 2 .6 kpc 
dShore et al.|[l990l ; Ide Winter & Thdfl99(il : lElia et alj|2004l) . Con- 
sequently its luminosity is not known either, being in the range 
1.9 < log(L/L ) < 4.9 (McGregor et al. 1988; Shore et al. 
1990, Hillenbrand et al. 1992). These uncertain stellar parameters 
and therefore the rather poor knowledge of the evolutionary phase 
of CD-42°11721 have stimulated us to investigate this object in 
more detail, based on a multiwavelength analysis. 
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Figure 1. The FEROS spectrum of CD-42° 11721. The most intense emission lines are indicated. A list of all identified lines is given in TablelATlin the 
Appendix. 



In this paper, we are presenting new high- and low-resolution 
optical spectra together with a more detailed description of avail- 
able infrared data taken from both the public catalogues as well 
as from the literature. We obtained a detailed optical spectral at- 
las, where the zoo of emission lines (many of them previously 
undetected) and the variety in line profiles give hints for the gas 
dynamics and distribution in the close-by circumstellar medium 
around CD-42°11721. We could also derive a set of stellar pa- 
rameters based on our high-resolution data, allowing to compare 
the position of CD-42° 1 1721 in the HR diagram with evolutionary 
tracks. Considering these stellar parameters, we have also modelled 
its SED considering different circumstellar scenarios: a spherically 
symmetric geometry, a passive flared disc and also an outflowing 
disc wind. 

The structure of the paper is as follows. In Sect.[2] we de- 
scribe our observations, the data reduction and the public data used. 
Sect.[3]gives a detailed optical spectral atlas and a description of the 



zoo of emission lines and the variety in line profiles found. In Sect. 
4 we present a description of the IR region. On the basis of our opti- 
cal spectra we derive in Sect.[5]the stellar parameters. In Sect. 6 we 
describe the SED modelling of this star, presenting the numerical 
codes developped by us, the results obtained and also the compari- 
son between them and the literature ones. In Sect. 7, we discuss the 
possible nature of CD-42° 1 1721 and present the conclusions of our 
study. 



2 OBSERVATIONS & REDUCTIONS 

Our optical spectra were obtained at the ESO 1.52-m telescope in 
La Silla (Chile) using the high resolution Fiber-fed Extended Range 
Optical Spectrograph (FEROS) and the low resolution Boiler & 
Chivens spectrograph. 

FEROS is a bench-mounted Echelle spectrograph with fibers, 
that cover a sky area of 2" of diameter, with a wavelength cov- 
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Figure 1 - continued 



erage from 3600 A to 9200 A and a mean spectral resolution of 
R = 48 000 corresponding to 2.2 pixels of 15 pirn. It has a com- 
plete automatic on-line reduction, which we adopted. Our FEROS 
spectrum presents the highest resolution used for observations of 
CD-42° 1 1721 to date. It was taken on June 10, 2000, with an ex- 
posure time of 5400 seconds. The S/N ratio in the 5500 A region 
is approximately 60. This high resolution spectrum and long wave- 
length coverage, shown in Fig.[T]will be used to obtain a detailed 
optical spectral atlas, for a detailed description of the different line 
profiles, to derive equivalent widths, and also to deblend some fea- 
tures in the low resolution spectra. 

The low resolution Boiler & Chivens spectrum (hereafter 
Cassegrain spectrum) was taken on the same observation run at 
ESO, but on June 9, 2000, with an exposure time of 900 sec- 
onds. The instrumental setup employed made use of grating #23 
with 600 1 mm - , providing a resolution of ~ 4.6 A in the range 
3800-8700 A. In the 5500 A continuum region, the S/N ratio of the 
Cassegrain spectrum is aproximately 180. Since there is no com- 



pletely line free region in the spectrum, the S/N derived is an upper 
limit. The slit width used has 2". Considering that the nebula in 
which CD-42° 1 1721 is embedded has dimensions of ~ 40" x 80" 
( Hutsemekers & Van Dromll 19901) . our observations were centered 
on the innermost parts of it. The Cassegrain spectrum was reduced 
using standard IRAF tasks, such as bias subtraction, flat-field nor- 
malization, and wavelength calibration. We have done absolute flux 
calibra tion using spectrophotometric standards from Hamuv et al.l 
Jl994l) . This spectrum is used to derive intensities of the emission 
lines for which we perform a detailed quantitative analysis. 

Equivalent widths and line intensities in the linearized spec- 
tra have been measured using the IRAF task that computes the line 
area above the adopted continuum. Uncertainties in our measure- 
ments come mainly from the position of the underlying continuum 
and we estimate the errors to be about 20 % for the weakest lines 
(equivalent widths lower than 2 A) and about 10 % for the strongest 
lines. 

We have also obtained photometric data from the litera- 
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ture, intending to describe the complete SED of CD-42° 11721. 
In addition, we have described the IR spectra obtained by 
ISO and IRAS satellites. These spectra are public and can be 
downloaded in the following web addresses: SWS01-ISO - 
http://isc.astro.cornell.edu/~sloan/library/swsatlas/aotl.html; 
LRS-IRAS - http://www.iras.ucalgary.ca/~volk/getlrs_plot.html 
The IRAS LRS spectrum covers a small wavelength range, from 
7.5 to 22 fim, with a spectral resolution of about 20-60. On the 
other hand, the ISO SWS01 spectrum has a higher resolution of ~ 
1500 and a wavelength coverage from 2 to 50 /im. 



3 THE OPTICAL SPECTRAL ATLAS: EMISSION LINES 
AND THEIR PROFILES 

To identify the lines and to create a spectral atlas of 
the optical region of CD-42° 1 1721, we have u s ed th e 
line lists provide d by iMoord dl945h; iThackeravl dl967h : 
lHamann & PerssorJ d 198a? : iMcKenna eta!] dl997l) and 



Landab errv. Pereira & de Arauiol feOOlh . We have also looked 
up two sites on the web: NIST Atomic Spectra Database Lines 
Form {http://physics.nist.gov/cgi-bin/AtData/lines_form) and The 
Atomic Line List v2.04 (http://-www.pa.uky. edu/^peter/atomic/\. 

FigQ] shows the complete high resolution spectrum of CD- 
42° 11721. It is completely dominated by emission lines. The only 
absorption lines visible are the Ca II and Na I lines, which are of 
interstellar origin, as well as the Diffuse Interstellar Bands (DIBs). 
Our spectra do not exhibit any photospheric lines. 

A huge number of permitted and forbidden emission lines 
were identified in the FEROS spectrum of CD-42° 1 1721 and listed 
in Table lAll given in the Appendix. The individual columns in Ta- 
ble lAll contain the observed wavelength (Col. 1), equivalent width 
(W(A), in A, Col. 2), the line profile type (Col. 3), and the pro- 
posed identification (Col. 4) for each line. The line identification 
given in Col. 4 encloses the ionization state of the element with the 
proposed transition and multiplet as well as the rest wavelength of 
the transition. It is possible that more than one ion can be allocated 
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to a single feature. In these cases, we give some possible alterna- 
tive identifications. For some lines no identification could be found. 
These lines remain unidentified, labelled as "Uid" in the tables. 

Inspection of Table lAl I reveals a huge zoo of emission lines 
from atoms mostly singly ionized, as well as a variety in their line 
profiles. Interestingly, there are also many emission lines from neu- 
tral metals like Cl, Nl, and Ol, and some sporadic traces for the 
emission of Si, and Mgl as well. Many of the lines listed in Ta- 
ble lAll are detected by the first time. 

In the following, we briefly discuss the presence of individual 
elements and lines and, if possible, compare our observations with 
previously published data. We especially highlight deviations of 
our observations from data in the literature which all seem to point 
in favour of either central source variability or ongoing changes in 
the CSM conditions. In addition, we classify the different line pro- 
files. 



3.1 The zoo of emission lines 

Hydrogen: The strongest lines in the spectra of CD-42° 11721 are 
the Balmer lines (Fig.[2}. We could identify the Balmer series up to 
H9 at the blue edge of our spectrum. While Ha, H/3 and Hy show 
indications for a double-peaked profile and broad wings extending 
to velocities of about 1400, 600, and 200 km s -1 respectively (for 
the Ha wings see Fig.[3}, the higher Balmer lines are singly peaked 
and their wing velocities decrease systematically. The equivalent 
width of Ha is almost a factor 2.5 lower than the value reported 
by Lopes, Damineli Neto & de Freitas Pacheco (19 92). The H/3 
equiv alent width has also been reported to be variable jLopes et al.l 
119921) . While it was found to be 27.7 A in the 1985 observations, it 
increased to 41.5 A in the 1988 observations. Our value, observed 
in 2000, lies with 30.4 A in between. 

In addition, at the long-wavelength end of the spectrum, the 
Paschen series arises. From this series we could clearly identify the 
lines from Pa(10) up to Pa(40). Their profiles are all narrow and 
single-peaked (see Fig.QJ. 

Helium: Helium seems to be variable in the spectrum of CD- 
42 1 17 21 . While it had (possibly) been detected in the years 1945- 
1948 bv | Merrill & Burwelll il949h , it was absent in the spectra of 
ICarlson & Henizd dl979h ~taken in 1949-1951 and 1962, and reap- 
peared again in the spectra observed in 1 985 bv lLopesetal]h992h 
and 1989 bv lHutsemekers & Van Droit] dl990f) . And also our data 
which have been taken in 2000 show emission of Hel. We do not 
detect any line from Hell. 

Oxyge n: The d e tectio n of the [Ol] 6300 A was first re- 
ported by lAcke et all d2005h - This line (and other oxygen lines) 
have not been seen durin g earlier observations perf ormed by 
Ide Winter & The] dl990h and lHutsemekers & Van Droml h990h . 

Our data , which have been taken two years earlier than those 
of lAckeet al. (2005), confirm the presence of this emission line, 
and we can report on numerous additional identifications of emis- 
sion lines of neutral oxygen, both forbidden and permitted. Inter- 
estingly, the [Ol] 6300 A and 6364 A lines show an asymmetric 
profile resembling a double-peaked line with a missing red peak, 
while the 5577 A line is clearly single-peaked and very narrow 
compared to the othe r two lines (Fig.[ H. The profile of the [Ol] 
line A 6300 shown bv lAcke et al.l J2005h agrees with ours, but there 
are significant dif ferences in the equivalent width of this line: while 
lAcke et ail d2005l) derived a value of 1.77 A, the equivalent width 
found from our observations is only 1.20 A. Since the set-up for 
both observations was quite similar, this difference in equivalent 
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Figure 2. Balmer lines from the CD-42° 11721 high resolution (FEROS) 
spectrum. The lines are broad and in emission with Ho, H/3 and Hy show- 
ing indications for a double-peaked profile. The absorption dip within the 
blue wing of He is the interstellar Ca II 3968 A line. Due to the scale of 
these figures, the extended wings of especially the Ho and H/3 lines cannot 
be seen. 
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Figure 3. The wings of the Ho line, extending to ~ 1400 km s~ 1 . The Cl 
A 6588 as well as the [N II] AA 6548, 6583 lines are also indicated. The 
dotted line represents the continuum. 



width might hint towards rapidly changing conditions, e.g. in den- 
sity and temperature structure, within the line forming region. 

In addition, we found (though very weak) the previously not 
detected [Oil] lines AA 7319, 7330. No single OlII line was identi- 
fied. 

Nitrogen: The only previously detected emission line from 
nitrogen in the spectrum of CD-4 2 1 17 21 is the line [Nil] A 5755, 
cited by Ide Winter & Thd h99fj) and iHutsemekers & Van Droir] 
(1990). The latter, however, ascribed the emission to the nebula 
and could not identify this line in the stellar spectrum. Our high- 
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Figure 4. The [O I] lines at 5577 A, 6300 A and 6363 A. Except the A 5577 
line, these lines are double peaked with a missing red peak. 

resolution spectrum confirms the existence of this emission line. In 
addition, we detected the lines [Nil] AA 6548, 6584, as well as Nil 
A 6173. Remarkably, we also identified many permitted emission 
lines from Nl. 

Carbon: The detection of any carbon line has not been re- 
ported yet. Our spectrum clearly shows the presence of several lines 
of Cl (see e.g. Fig.[3}, one line of ClI and Cm each, as well as the 
possible detection of a [Cl] line. 

Sulfur: The existe nce of the [Sll] AA 6716, 6731 line s has for- 
merly been reported by Hutsemekers & Van Droml dl99fj|) as com- 
ing from the nebula. In addition to these lines, our spectrum also 
contains the line [Sll] A 4068 while the A 4076 line is missing. 
From Sill we could identify only one line, i.e. [Sill] A 9069. The 
line [Sill] A 631 1 does not show up in our spectrum. However, as in 
the case of N, O, and C, we have some indications for the presence 
of Si emission. 

Iron: It is the element with by far the largest number of emi- 
sison lines. With the exception of four lines allocated to Fein (two 
forbidden and two permitted), all lines are from singly ionized iron. 
CD-42 11721 is known t o show a huge number of Fell lin es (e.g. 
ICarlson&Henizelll979l : Hutsemekers & Van Droml 1 19901) . how- 
ever, the existence of an equally large amount of [Fell] lines as we 
found in our spectrum, has previously not been reported. Actually 
Ide Winter & The! J 199Ct) reported on a possible detection of a [Fell] 
line, however this identification is rather doubtful due to their low 
spectral resolution. These lines show either single-peaked, double- 
peaked, or even multiple-peaked profiles, indicating a high com- 
plexity of the CSM. The permitted lines are usually more intense 
than the forbidden ones and also broader with wings extending to 
100-200 km s"\ while the line wings of the forbidden lines indi- 
cate velocities of about ~ 50 — 60kms _1 , only. 

Other elements: In addition, we identified emission lines 
from the following metal ions: Cm, Mgl, Mgll, Mnll, Till, Sill, 
and Nell. Most of them have previously not been detected. 

3.2 Classification of the line profiles 

To understand the structure of the CSM of CD-42° 11721, we sort 
the observed emission line profiles into three categories (see Ta- 
ble lAlb : (i) single-peaked, (ii) double-peaked, and (iii) multiple- 
peaked. In agreem ent with the descriptions in earlier studies (e.g. 
lLopes etal.ll 19921) . no single P Cygni profile was detected. 

Single-peaked profiles: These profiles clearly dominate. 
They can be found without any exception for every element in ev- 



Table 1. Velocities of the H Balmer lines and some forbidden lines seen in 
the FEROS spectrum. Listed are the ions, the laboratory wavelength of the 
line, the blue and red wing velocities, the mean expansion velocity, and the 
line center velocity. 
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Table 2. Peak separation of double-peaked lines for different ions. 
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ery ionization stage, independent whether these are forbidden or 
permitted lines. To draw any conclusion about the dynamics of the 
CSM we measured the line wing velocities of several forbidden 
(and therefore optically thin) emission lines. The results are listed 
in TableQ] together with the high wing velocities found from the 
Balmer lines. Two striking conclusions can be drawn from this ta- 
ble: (i) the mean wing velocity (which we call the expansion veloc- 
ity, v exp ) is on the order of 50-60 km s _1 and about constant for 
all the forbidden lines, and (ii) the line center velocity of all lines is 
also about constant and on the order of -20 kms" 1 . 

Double-peaked profiles: There are quite a lot of lines show- 
ing a double-peaked profile. The blue and red peaks are thereby of 
equal strength with only some rare exceptions. One remarkable ex- 
ception are the [Ol] lines at 6300 A and 6364 A for which the red 
peak is missing (Fig.[4]l. Another exception are the Balmer lines, 
especially Ha and H/3 (see Fig.[2jl 

In Table|2]we list some of the well resolved double-peaked 
lines from the FEROS spectrum together with their peak separation, 
which lies between 30 and 50 km s" 1 (neglecting the Balmer lines). 

Interestingly, the double-peaked profiles are present only for 
neutral metals (like Ol, Nl, Cl, and Si) and ions with low ionization 
potential, i.e. lines from Fell, Mgll, and Cm with \ < 8 eV, which 
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is well below the ionization potential of hydrogen. This fact has 
two severe implications: 

• the emitting material must be neutral in hydrogen, and 

• the double-peaked profiles might indicate rotation. 

Combination of both seems to suggest the existence of a neu- 
tral (in hydrogen) rotating disc or an equatorial outflowing wind 
(see e.g. Zickgraf 2003). However, we emphasize that the double- 
peaked profiles are not restricted to permitted lines only, but are 
also present in several forbidden lines from [Ol], [Cm], and of 
course [Fell]. 

Multiple-peaked profiles: These types of profiles are found 
especially for permitted emission lines of Fell where three peaks 
or more are observed. The existence of such multiple-peaked lines 
favours a highly complex multi-component emission region. 

3.3 Comparison with optical spectra of other B[e] stars 

Based on low-resolution optical observations, it has been noted in 
the literature that the spectrum of CD-42° 1 1721 is quite similar to 
some other B[e] stars. One of these is HD 87643 (Hutsemekers & 
van Drom 1990; Oudmaijer et al. 1998), another object belonging 
to the group of unclassified B[e] stars since it also shows charac- 
teristics that might cla ssify this object e ither as a HAeB[e] stars or 
as a B[e] supergiant dLamers et al. 1998). The low-resolution spec- 
tra as shown in Fig.[5]are indeed similar, however, there exist also 
differences that show up already in this plot: in contrast to CD- 
42 11721, the spectrum of HD 87643 clearly shows the presence 
of P Cygni profiles. 

Another object found to exhibit a similar optical low- 
resolution spectrum (also included in Fig.[5j, is CPD-52°9243 
(Carlson & Henize 1979; Swings 1 981; Winkler & Wo lf 1989), a 
galactic B[e] supergiant candidate fcamers et alj|l998h . The most 
obvious difference in its optical spectrum is, however, the presence 
of the Hel lines in absorption. 

For these three stars, high-resolution (FEROS) observations 
show very clearly that there exist numerous other differences (see 
Borges Fernandes 2004), which cannot be listed here in detail. 



4 THE INFRARED REGION: GAS AND DUST 
FEATURES 

Concerning to the IR region, CD-42° 11721 was observed by ISO 
and IRAS satellites. The IRAS LRS spectrum (Fig. 6) is character- 
ized by a strong absorption centered around 10 /im, originated by 
amorphous silicate, that is not seen in the ISO SWS01 spectra. It 
was classified by Zhang, Chen & He (2004) as "noisy or incom- 
plete", meaning that the absorption could not be real. Actually, the 
first part of LRS IRAS (from 7.6 to 13.4 /im) is not so noisy and we 
can clearly see the previously known as "Unindentified IR bands" 
(UIR) that are commonly accepted nowadays as being caused by vi- 
brational transitions of Polycyclic Aromatic Hydrocarbons (PAH), 
at 7.6, 11.3 and 12.7 pm. 

The SWS01 ISO spectrum (see also Fig. 6) of this object 
shows especially a strong rising continuum beyond 15 \im, being 
clearly seen gas and dust emission features (Voors 1999). Concern- 
ing the gas emission is remarkable the presence of many H recom- 
bination lines (Benedettini et al. 1998) and also the presence of 
forbidden lines, e.g., [Sin], [Fell], and [S III] lines. Other forbid- 
den lines, e.g., [Ol], [Cll], [Nil], and [Oin] lines, were identified 
by Lorenzetti et al. (1999) by the analysis of the Long Wavelength 
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Figure 5. Comparison of the Cassegrain spectra of CD-42° 11721 with 
some other B[e] stars. Though quite similar at first glance, they show severe 
differences: the spectrum of HD 87643 does show the presence of P Cygni 
profiles, and the spectrum of CPD-52° 9243 has the Hel lines in absorption. 

Spectrograph (LWS-ISO) spectra of this star. We have also tenta- 
tively identified some features as being caused by H2 molecule, 
however these identifications are very uncertain. McGregor et al. 
( 1988) cited the probable presence of CO bands, however we could 
not identify these features in the SWS01 ISO spectrum. 

Regarding the dust emission, the ISO-SWS spectrum presents 
the called "dual-dust" chemistry, characterized by the existence of 
solid state bands of C-rich and O-rich dust in the same environ- 
ment. The C-rich material is represented by PAH emission bands 
(lourdain de Muizon, d'Hendecourt & Geballe 1990; Voors 1999). 
On the other hand, beyond 20 [im, the O-rich dust is represented 
by some features caused by crystalline silicates. 

In Table 3, we compile the identifications available from dif- 
ferent papers and those ones determined by this work, giving a 
complete view of the most important features present in the ISO- 
SWS spectrum of CD-42° 1 1721, from 2 to 50 pm. There, the ob- 
served wavelength (Col. 1), the proposed identification (Col. 2), 
and the reference (Col. 3) for each feature is given. 

However, all the spectral description cited above, in this sec- 
tion, probably is not related to the star itself. This happens be- 
cause ISO and IRAS data cover a large area of the sky around 
CD-42° 1 1721 and therefore they are probably contaminated by the 
known surrounding reflection nebula and/or other stars, since it lies 
in a crowded region (Habart et al. 2003; Hamaguchi et al. 2005). 
This contamination is clearly shown by Acke & van den Ancker 
(2006) that presented new IR spectra of the 3 and 10 pm wave- 
length regions. These spectra cover an area of only 1" of diameter 
around CD-42° 11721, therefore the region close-by this star, and 
surprisingly, they do not present the strong PAH emission seen in 
the ISO and IRAS spectra. This fact means that "dual-dust" chem- 
istry is not related to the object itself. 
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Figure 6. The SWS01 ISO spectrum (solid line) and the IRAS LRS spec- 
tram (dotted line) of CD-42° 1 1721 . 

5 THE STELLAR PARAMETERS OF CD-42° 11721 

Most of the problems in classifying CD-42° 11721 are caused by 
the uncertain stellar parameters. We therefore use firstly our optical 
high-resolution data for a tentative determination of the stellar pa- 
rameters. The details on how we estimate each parameter are given 
in this section. 



5.1 The effective temperature 

The effective temperature values f or CD-42° 11721 found in the 
literature vary between 12 300 K [Hille nbrand et alj fl992l) and 
31600K (Cidale et al. 2001). Since in our spectra there is the 
presence of emission Hel lines and the absence of Hell lines, we 
can therefore conclude that the effective temperature must be in 
the range between 13 000 K (see Machado & de Araujo 2003) and 
30 000 K. 

IUE spectra have been used to constrain the spectral type of 
CD-42° 1 1721 to early B (or B0) by Shore et al. (1990). This result 
should however be taken with caution, because: (i) CD-42° 11721 
lies in a crowded region with even some close-by X-ray sources 
(see Hamaguchi et al. 2005), (ii) with an elliptical aperture with ex- 
tensions of 10" x 20", IUE observed definitely more than one star, 
and (iii) the IUE flux is much stronger and clearly offset from our 
optical spectrum, also indicating that multiple sources have been 
observed. The IUE spectrum can therefore not be used for a proper 
temperature determination and the classification as B0 (or early B- 
type) star is therefore not justified. 

The absence or presence of individual ionization stages of the 
metals identified in our spectrum (see Table lAlt can be also used 
for a better restriction of the effective temperature. Inspecting the 
ionization potentials for the different ions with observed lines in 
our spectra, we can state that CD-42° 11721 shows clear hints for 
the presence of ions with an ionization potential below \ — 25 eV 
like e.g. Cm, Oil, Nil, Nell, Hel recombination lines, and (even 
though very weak) Sill. On the other hand, no lines from ions with 
ionization potential higher than ~ 27 eV could be identified. Good 
tracers for a higher ionization potential are the lines from Ann, 
NlII, and Clin as have been obs erved e.g. from th e compact plan- 
etary nebula B[e] star Hen 2-90 dKraus et"al]|2005h . From this nar- 
row range between 25 eV and 27 eV, which is just above the Hel 
ionization potential, we conclude that the effective temperature can 
only slightly exceed the value of 13 000 K (Machado & de Araujo 



Table 3. The identification of the features present in the ISO-SWS spec- 
trum. There Ident. means the proposed feature identifications, enst. is en- 
statite, diop. is diopside, fost. is fosterite, meth. gr. is related to me- 
thil groups, and Uid means unidentified features. References (Ref.): B = 
Benedittini et al. (1998), V = Voors (1999). It is important to cite that ex- 
cept where noted, otherwise, the identifications are ours. 
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2003). In fact, in order to reproduce the shape of the optical contin- 
uum (which for a given extinction is strongly temperature sensitive, 
see Sect. 5.4), the effective temperature cannot exceed the value of 
about 15 000 K. We thus conclude that CD-42° 1 1721 has an effec- 
tive temperature of 14 000 K ± 1 000 K. 



5.2 The distance 

The most controversal parameter of CD-42 1 1721 is its distance. 
Literature values range from 136pc dElia et alj|2004l) over 160 pc 
dHillenbrand et al J 19921) and 400 pc dde Winter & The]|l990t) up to 
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Table 4. Observed line ratios used for the extinction derivation. 



n 


Pa(n)/He 


Pa(n)/H<5 


Pa(n)/H7 


17 


1.308 


0.662 


0.357 


20 


0.876 


0.443 


0.239 


21 


0.799 


0.404 


0.218 


24 


0.476 


0.241 


0.130 



2 kpc ( McGregor et al. 19 88), 2.5 kpc dLopes etal]|l992h and even 
2.6 kpc ( Sho re et aljl990h . This range in distances by about a factor 
20 leads to a scatter in stellar luminosity by a factor 400. 

For the distance determination from our spectra we make use 
of the relati ons for the i nterstellar absorption lines of Nal and Call- 
K given by Allen] a 19730 - From the equivalent widths of these lines 
(see Table lAlt we find the following distances: 

• d = 1.23 kpc derived from the Nal lines, and 

• d = 1.06 kpc derived from the Call-K lines 

Both values agree quite well and we will adopt the mean value 
of d ~ 1.15 kpc, which lies between the very small and very large 
distances found in the literature, as a reasonable value. Based on the 
relations used, we can estimate the error in distance. Taking into 
account the errors in the measurements of the equivalent widths, 
on the order of 20 %, the error in distance results in 0.15 kpc. We 
are aware that this error does not reflect the full error, that might 
in principle be larger than this, since the distance estimated from 
the use of interestellar absorption lines contains a systematic error 
due to the statistical relations used. In our opinion, this parameter 
remains as the worst estimated, consequently all the conclusions 
that depend directly on it must be taken with caution. 

5.3 The interstellar extinction 

The extinction towards CD-42° 1 1721 is believed to be quite high. 
Hints for this are provided by the strength of the interst ellar diffuse 
featur e at 4430 A and the faint IUE SWP spectrum ( Sh ore et all 
1990). Literature values for the visual interstellar extinction, Av, 
are found to range from 4.2 mag (Cidale et al. 2001) to 7.1 mag 
dde Winter & Thelll990l) . With such an enormous range of extinc- 
tion values, it is impossible to draw any reliable conclusion from 
our optical spectra. We therefore attempt to further constrain the 

extinction. 

As a first attempt we tried to use the Pagel method jPagell 

1969). This method is based on a relation between the observed 
fluxes of Fe II forbidden lines and the reciprocal wavelengths. How- 
ever, in the case of CD-42° 1 1721, it could not be employed due to 
the lack of enough non-blended [Fe II] lines in our high-resolution 

spectra. 

W e therefore make use of the method suggested by iKrauj 
d2007l) that is based on the line ratios of the Paschen lines over the 
Balmer lines. The red edge of our FEROS spectrum shows the pres- 
ence of well resolved Paschen lines from Pa(10) up to Pa(40). In 
addition, on the blue edge of the spectrum, we observe the Balmer 
lines up to H9. T he most unaffe cted line ratio to derive the extinc- 
tion according to iKrausI d2007t) is the line ratio Pa(10)/H9. Since 
the FEROS data are not flux calibrated, we derive the line strengths 
by calibrating the data with the Cassegrain spectrum. This spec- 
trum has, however, a narrower wavelength range than FEROS. We 
extrapolated the continuum of the Cassegrain spectrum to the blue 
and red edges of FEROS, but the lines in these ranges have a much 



Table 5. Mean extinction values, Ay in magnitudes, from the line ratios 
calculated for different mass fluxes and electron temperatures. The Paschen 
lines used to calculate the mean values are Pa(17), Pa(20), Pa(21), and 
Pa(24). 
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Pa(n)/H5 


12 000 
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4.571 


4.298 
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Pa(n)/H7 


12 000 


5.153 


4.688 


4.493 


4.521 


Pa(n)/He 


15 000 
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4.697 


4.409 


4.344 


Pa(ra)/H<5 


15 000 


4.855 


4.464 


4.208 


4.197 


Pa(n)/H7 


15 000 


4.984 


4.535 


4.354 


4.393 


< A v > 


12 000 


5.059 


4.678 


4.422 


4.398 


< A v > 


15 000 


4.945 


4.565 


4.324 


4.311 



larger uncertainty and are therefore not so useful to derive the 
reddening. We therefore concentrate on the line ratios Pa(n)/He, 
Pa(n)/H<5, Pa(n)/H7 with n = 15 . . . 25. The reasons for these re- 
strictions are the following: 

• Paschen lines from levels higher than n = 25 are weak and a 
proper determination of their line strengths becomes difficult. 

• Paschen lines from levels lower than n = 15 suffer from 
the uncertain continuum extrapolation beyond the red edge of the 
Cassegrain spectrum. 

• The H9 line is affected by the uncertain continuum extrapola- 
tion of the Cassegrain spectrum which has a rather noisy continuum 
at the blue edge. 

• The H8 line is not useful since it is blended with a Hel line. 

Restricting to unblended, clearly resolved Paschen lines (see 
Fig.QJ reduces the lines suitable for a proper extinction determi- 
nation to Pa(17), Pa(20), Pa(21), and Pa(24). For these lines the 
observed line ratios with the Balmer lines He, H5, and H7 are sum- 
marized in Table|4] 

To determine the extinction, we need to know (at least approx- 
imately) the electron temperature in the wind and the mass flux of 
the star. With an effective stellar temperature of about 15 000 K and 
the fact that in a stellar wind the Paschen lines form very close to 
the central star where the wind temperature is still rather high, we 
use the unreddened line ratios for the two different electron tem- 
peratures, 12000K and 15 000 K. The mass flux of the star is not 
known a priori. But since the different line ratios should result in 
the same (or a very similar) extinction value, we can exclude the 
ve ry high mass fluxes (i.e. F m ^ 10~ 4 g s _1 cm -2 ) in the tables 
of lKrausI d2007f) , because they result in a large spread in extinction 
values. 

To pin down the extinction, we first calculate a mean value 
of the different line ratios found for each mass flux in the range 
Fjn — 10 -5 . . . 10~ 8 g s _1 cm~ 2 and for the two different electron 
temperatures. These mean values are given in Table[5] From this 
table it is clear already, that the visual extinction must be in the 
range from 4.3 mag to 5.0 mag, which is in rough agreement with 
the values of 4.2 mag and 4 .3 mag found by Cidale et al. (2001) 
and lHillenbrand et all dl992h . and the value of about 5 mag given 
by McGregor et al. (1988). 

To restrict the extinction even further, we next calculate the re- 
combination line luminosities of the Paschen lines, using the model 
parameters for T e and F m , redden them with the corresponding ex- 
tinction value found, and compare them with the observed line lu- 
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Table 6. Stellar parameters for CD-42° 11721 found from our high- 
resolution optical spectra as described in Sect.[5] The mass and log g values 
have been derived under the assumption that CD-42° 11721 is a post-main 
sequence object. The values in the second row represent the error in each 
measurement. 



T aH R, L* M initial log g d A v 

[K] [R ] [Lq] [M ] [kpc] [mag] 

14000 17.3 1.0 X 10 4 8-10 ~3 1.15 4.8 

±1000 ±0.6 ±0.3 x 10 4 ±0.15 ±0.2 



5000 6000 7000 

Wavelength (A) 

Figure 7. Fitting of the observed op tical co ntinuum (dashed line) with Ku- 
rucz model atmospheres (solid lines, Kurucz 1979). Shown is one of the best 
fits achieved for a solar metalicity star with T e ff = 14 000 K, log g = 3.0, 
and a foreground extinction of Ay = 4.8 mag. 



minosities under the assumption of a distance of 1.15 kpc. Such a 
comparison can give only a rough result, because the line lumi- 
nosities, different from the line ratios, do depend on individual 
stellar parameters like the stellar radius and mass loss rate. As- 
suming spherical symmetry, which should result in a lower limit 
of the real mass flux of the star, and varying the stellar parameters 
in a reasonable range, we found that even with this rough method 
we can exclude the mass fluxes lower than 10~ 6 gs~ 4 cm" 2 . This 
means, that the extinction towards CD-42° 1 1721 is on the order of 
4.8 ±0.2 mag. 
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Figure 8. The SED of CD-42° 11721, considering all data with different 
apertures presented in Table 7. 



5.4 Stellar radius and luminosity 

The stellar radius of CD-42° 11721 can be found by fitting the 
Cassegrain sp ectrum with classical Kurucz model atmospheres 
jKurucz|[l979t) . Since the effective temperature of the star is in the 
range 13 000 to 15 000 K, we fitted Kurucz model atmospheres for 
stars with solar metallicity and different log g values, reddened with 
4.8 ± 0.2 mag interstellar extinction. From these fittings it turned 
out that no reliable fit was found for Av = 5.0 mag, while for 
the values of Av — 4.6 — 4.8, especially the models with low to 
mid surface gravity, fitted best (Fig. 7). From the fitting parame- 
ter which is simply R^/d 2 and for a distance of 1.15 ± 0.15kpc 
we can derive the stellar radius. This has been done for all the best 
fit model atmospheres. The resulting stellar radius taking into ac- 
count the error in distance is i?» = 17.3 ± 0.6 Rq. We want to 
stress that this is an upper limit for the stellar radius, because the 
optical continuum might contain some contributions from free-free 
and free-bound emission generated in the wind. Nevertheless, the 
contribution of the wind can only be of minor impact because the 
photometric data in the near infrared do not show influences of 
free-free emission (see Sect. 6.1). We therefore regard the derived 
~ 17 R© as a reasonable value for the stellar radius. Since we know 
the radius and the effective temperature, we could find a stellar lu- 
minosity of L, = (1.0 ±0.3) x 1O 4 L . 

To finish our optical analysis, we summarize in Table 6 the 
stellar parameters of CD-42° 1 1721 derived within this section. 



6 MODELLING OF THE SED OF CD-42° 11721 

To construct a proper SED for CD-42° 1 1721, we collected all pho- 
tometric data available in the literature and listed them in Table 7. 
There, the observed wavelengths (Col. 1), the measured fluxes (Col. 
2), the aperture diameters (Col. 3), and the references (Col. 4) are 
given. As can be seen in Fig. 8, considering all data and different 
apertures, the SED is double peaked, similar to those HAeBe stars 
called "group I" (Meeus et al. 2001), though differently than other 
objects of this group, its first peak is shifted to higher wavelengths, 
around 2 p_m. However, as cited in Sect. 4, the data from observa- 
tions with large apertures are contaminated by the reflection nebula 
and other sources and do not come from the star itself. Based on 
this, we have decided to consider in our modelling only data with 
aperture diameters equal or smaller than 15". It is important to cite 
that the most confident results would be obtained from the mod- 
elling of a SED considering only photometric data using apertures 
around 1" or 2", however there is not enough available data to build 
a reliable SED for this object. 

In this section, we present the modelling of the SED of CD- 
42° 11721, using three diferent codes, that solve the problem of 
radiative transfer in a circumstellar dust medium, considering three 
different scenarios, e.g. a spherically symmetric geometry, a pas- 
sive flared disc and an outflowing disc wind. We aim to determine 
which scenario describes better the SED of this star, assuming the 
set of stellar parameters obtained from our optical analysis. Our 
codes that represent the flared disc and the outflowing disc wind 
scenario are presented here for the first time. 



Table 7. Photometric data of CD-42° 11721 available in the literature. 
Refs: 1 = de Winter et al. (2001); 2 = Herbst (1975); 3 = Berrilli et al. 
(1992); 4 = Elia et al. (2005); 5 = Mannings (1994); 6 = Hillenbrand et al. 
(1992); 7 = Henning et al. (1994); 8 = Skinner, Brown & Stewart (1993); 
ISO-SWS = Acke & van den Ancker (2004); MSX = Egan et al. (1999); 
IRAS = http://irsa.ipac.caltech.edu/IRASdocs/iras_mission.html , 2MASS = 
http://www.ipac.caltech.edu/2mass/releases/allsky/doc/secl_6b.html 



Wavelength Flux Aperture Diameter* References 

[|_m] [iym -2 |_m _1 ] [arcsec] 
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Table 7 - continued 



Wavelength 


Flux 


Aperture Diameter* References 


[pin] 


[VKm _2 |_m _1 ] 


[arcsec] 


1300.00 


7.54xl0~ 19 


23.0 7 


36000.00 


3.06xl0" 24 


2.5 8 



* = except by IRAS and ISO-SWS, that have rectangular apertures. ** = 
Mean values of observations with same aperture and from the same refer- 
ence. 



6.1 Free-free and free-bound emission 

For the modeling of the SED of CD-42° 1 1721, our first step was to 
determine the contribution of the free-free and free-bound emission 
(hereafter ff-fb emission) to the total continuum. The ff-fb con- 
tinuum can be calculated by accounting for the observed 3.6 cm 
emission that was obtained with an aperture size of 2.5" (see Ta- 
ble 7) and that can be ascribed purely to the ionized wind ma- 
terial. Since we do not know exactly the wind geometry and the 
possible inclination under which the system of CD-42° 11721 is 
seen, we calculated the ff-fb continuum assuming an isothermal 
(T e = lOOOOiO spherically symmetric wind, with an outer edge 
of 1500 AU. The mass loss rate needed to fit the 3.6 cm emission 
(Al — 2.86 x 10~ 6 MQt/r _1 ) must be considered as a lower limit. 
The terminal velocity, which is proportional to the escape velocity 
on the stellar surface, was calculated from the stellar parameters 
given in Table 6 and it is 515fcms _1 . It turns out, that the contribu- 
tion of the ff-fb continuum is negligible in the optical and near-IR 
spectral range (see Fig. 9). 



6.2 Spherically symmetric envelope 

Our optical results revealed that the circumstellar scenario of CD- 
42° 11721 should probably be composed by a non spherical geom- 
etry. However, for sake of completeness, we modeled the SED of 
this star considering also a circumstellar spherically symmetric dust 
envelope. 

Our numerical treatment of radiative transfer in a spherical 
envelope was described in Lorenz-Martins & Pompeia (2000). The 
Monte Carlo scheme is used for representing the propagation of ra- 
diative energy photon by photon. For each interaction between a 
"photon" and a grain, a fraction of the energy is stored (absorption) 
and the remaining part is scattered according to the scattering di- 
agram. Two species of grains are used and they interact with each 
other. The stellar radiation leads to a first distribution of dust tem- 
perature and the thermal radiation from grains is simulated, giv- 
ing after several iterations the equilibrium temperature. Here two 
species of grains are considered both with a numerical density de- 
creasing as r~ 2 , which corresponds to an expansion at constant ve- 
locity. The grains can be present at same or at different distances of 
the central star. The following physical quantities are required in a 
first guess and are subsequently fixed by fitting the infrared/near in- 
frared flux: the effective temperature of the central star, T e ff\ grain 
sizes and refractive indices for all wavelengths; the extinction opac- 
ity, T ex t at A = 9 fim; the inner radii of the dust envelope (say, Rsu, 
inner radius for amorphous silicate grains; Rcrys, inner radius for 
crystalline silicate grains, or Ra.c. amorphous carbon grains), and 
R the outer radius of the dust envelope. 

The models give not only the spectral repartition of the to- 
tal flux but also of its different components (direct, scattered, ther- 
mal for each species of grain), as well as the angular distribution 



0.36 


4.67xl0~ 13 


15.0 


1** 


0.36 


3.86xl0~ 13 


9-13 


2 


0.44 


6.23xl0~ 13 


15.0 


1 ** 


0.44 


5.44xl0~ 13 


9-13 


2 


0.56 


1.16xl0 -12 


21.5 




0.56 


1.07xl0~ 12 


15.0 


j * * 


0.56 


9.22xl0~ 13 


9-13 


2 


0.64 


2.17xl0~ 12 


15.0 


1 ** 


0.79 


2.66xl0~ 12 


15.0 


1** 


1.23 


3.69xl0 -12 


15.0 


3 


1.25 


4.40xl0~ 12 


13.0 


1*» 


1.25 


3.92xl0 -12 


8.0 


2MASS 


1.63 


5.11xl0 -12 


15.0 


3 


1.65 


5.37xl0~ 12 


13.0 


1*» 


1.65 


4.81 xlO -12 


8.0 


2MASS 


2.19 


6 60x 10 -12 


15.0 


3 


2.20 


7.11 x 10 -12 


13.0 




2.20 


6.33xl0~ 12 


8.0 


2MASS 


3.30 


6.50xl0~ 12 


14 x 20 


ISO-SWS 


3.40 


6.49x 10 -12 


14 x 20 


ISO-SWS 


3.45 


8.31xl0 -12 


13.0 




3.50 


6.29x 10~ 12 


14 x 20 


ISO-SWS 


3.79 


5.83x 10~ 12 


15.0 


3 


4.29 


6.86xl0~ 12 


18.3 


MSX 


4.35 


6.89x 10 -12 


18.3 


MSX 


4.64 


4.87x 10~ 12 


15.0 


3 


4.80 


4.37x 10~ 12 


13.0 


1** 


6.20 


3.23xl0~ 12 


14 x 20 


ISO-SWS 


7.70 


2.24xl0~ 12 


14 x 20 


ISO-SWS 


8.28 


3.26xl0~ 12 


18.3 


MSX 


8.38 


1.73xl0~ 12 


15.0 


3 


8.60 


1.77x 10~ 12 


14 x 20 


ISO-SWS 


9.69 


9.06 x 10 -13 


15.0 


3 


9.70 


1.26x 10~ 12 


14 x 20 


ISO-SWS 


1 1 00 


1 14x 10 — 12 


14 x 20 


low .3 vv . > 


12.00 


1.98xl0~ 12 


45.0x270.0 


IRAS 


12.00 


9.26x10" 13 


80.0 


4 


12.13 


1.54xl0~ 12 


18.3 


MSX 


12.89 


4.90xl0~ 13 


15.0 


3 


14.65 


6.47xl0~ 13 


18.3 


MSX 


21.34 


7.84xl0~ 13 


18.3 


MSX 


25.00 


1.21X10" 12 


45.0x276.0 


IRAS 


25.00 


6.68x10" 13 


80.0 


4 


60.00 


1.60xl0~ 12 


90.0x282.0 


IRAS 


60.00 


1.31xl0~ 12 


80.0 


4 


100.00 


6.63x10" 13 


180.0x300.0 


IRAS 


100.00 


5.12xl0~ 13 


80.0 


4 


350.00 


2.32X10- 16 


18.5 


5 


450.00 


8.02xl0~ 17 


18.3 


5 


600.00 


l.llxlO- 17 


17.5 


5 


750.00 


4.74x10" 18 


18.0 


5 


800.00 


5.06xl0~ 18 


17.0 


5 


850.00 


2.86xl0" 18 


18.0 


5 


1 100.00 


6.20xl0~ 19 


18.8 


5 


1300.00 


4.58x10" 19 


19.8 


5 


1300.00 


4.79xl0~ 19 


28.0 


6 



12 Borges Fernandes, Kraus, Lorenz Martins & de Araujo 



Table 8. Results of the modelling considering a spherical symmetric enve- 
lope. The temperatures are given in K. More details are in the text. 



Model 


J eff 


Ta.c. 


Tsil 




Tcrys 2 


A 


10000 


564 


364 






B 


15000 


781 


570 






C 


25000 


1038 


880 






D 


10000 






346 


327 



of these radiation fields and the radial dependence of dust tem- 
perature. The absorption and scattering efficiencies, as well as the 
albedo for the grains were calculated by us using the Mie theory 
and the optical constants given by Rouleau & Martin (1991), David 
& Pegourie (1995) and Jaeger et al. (1998). 

Several models, considering different sets of parameters, in- 
cluding those derived by our optical analysis, were calculated by 
us. The results of four of them are presented in Table 8. There, the 
name of each model (Col. 1) is followed by T e // (in K, Col. 2) 
and the temperatures (also in K) of the hottest grains (Cols. 3, 4, 5 
and 6). Other parameters are kept constant in these models. These 
are r ext = 3.0 at A = 9, Rsu = Rcrys = Ra.c. = 2007?* and R c 
= 5000i?* . The grains sizes are also fixed at 200A and the ratio of 
amorphous carbon over amorphous silicate grains is set to 0.50 and 
the ratio of crystalline silicate grains needed in model D is also set 
to 0.5. 

As can be seen in Fig. 9, the total SEDs (including the stellar 
and dust contributions) obtained by us, assuming a spherical geom- 
etry for the dust envelope, could not model properly the SED of 
CD-42° 11721, especially the optical region. The model A and D 
differ one from the other only by the dust species (see Table 8). In 
model A we have used amorphous silicates and amorphous carbon 
grains simultaneously, while in model D we have considered only 
crystalline silicates. The other two models (B and C) differ from 
the model A only by T e f / as can be seen in Table 8. These models 
present very similar fitting, even considering quite different set of 
parameters. 



6.3 Equatorial disc 

Based on our optical analysis, a more complex circumstellar geom- 
etry, most plausibly a disc scenario, seems to be indeed necessary 
to account for the observed SED of CD-42° 1 1721. 

The presence of a disc seems to be the better answer to explain 
the spectral characteristics of not only B[e] supergiants (Zickgraf 
et al. 1985), but also of HAeBe stars. The discs around young stars 
are detected with millimeter interferometry (Mannings & Sargent 
1997, 2000; Fuente et al. 2004), and in some cases, they are visible 
in optical images of the stars, either via scattered light or in silhou- 
ette against a background light source (Burrows et al. 1996; Grady 
et al. 2000; McCaughrean & O'Dell 1996). 

Due to the absence of images and also interferometric obser- 
vations for CD-42° 11721, we have decided to reproduce its SED 
considering two different types of discs, one typical of HAeBe stars 
and another one, typical of B[e]sg stars. 




10 100 1000 in" 10 5 

Wavelength (jum) 

Figure 9. Different models reproducing the SED of CD-42° 11721, consid- 
ering a spherical symmetric envelope. We can see the model A (dashed 
line), model B (long-dashed line), model C (dot-short dashed line) and 
model D (dotted line). A detailed explanation of these models can be seen 
in the text and also in Table 8. The triangles represent the photometric data 
obtained in the literature. The dot-long-dashed line is the ff-fb emission that 
fits the 3.6 cm data point. Its contribution, as can be seen, to the IR part of 
the SED is negligible. 



6.3.1 Passive flared disc 

Considering CD-42° 11721 as an HAeBe star, we need to assume 
the existence of a passive equatorial disc, where passive means that 
the disc is heated purely by stellar irradiation. 

Specifically for CD-42° 11721, we need to explain the loca- 
tion of the hot dust that might be responsible for the peak seen in 
its SED. We are assuming here that the hot dust is located at the 
surface of the circumstellar disc. For the surface to become much 
hotter than the interior, it must be inclined with respect to the in- 
falling stellar radiation, i.e. the disc must be flared. Such a flaring 
can happen naturally when the disc is in hydrostatic equilibrium in 
vertical direction with a radial temperature distribution that falls off 
more slowly than r _1 (Kenyon & Hartmann 1987). 

The flared disc is proposed to explain the spectral character- 
istics of young pre-main sequence stars of the "group I" (Meeus 
et al. 2001; Dullemond & Dominik 2004; Acke et al. 2005). The 
geometry of this kind of disc also explains the presence of a more 
intense flux in the far-IR, seen in HAeBe stars of the "group I", 
since the most distant regions in the disc surface receive a higher 
amount of stellar radiation (they see a larger solid angle of the star). 
This kind of disc can also explain the presence of the "dual-dust" 
chemistry, since the PAH molecule emission would be formed in 
the surface layers of the disc and the silicate bands in the inner 
regions of it, closer to the star (Meeus et al. 2001; Dullemond & 
Dominik 2004). However, as cited in the Sect. 4, in the case of CD- 
42° 11721, this mixed chemistry probably does not belong to the 
circumstellar medium of this star. 

Considering a flared disc scenario, we developped an analyt- 
ical two-layer model that contains an optically thick (r™ ld 3> 1) 
mid-layer which is isothermal in ^-direction. The mid-layer is sand- 
wiched by two top-layers which are marginally optically thick at 
visual wavelengths but still optically thin at IR wavelengths, i.e. 
t"v° P ^ 1 and r^ p ^ 1. Consequently, the infalling stellar 
light is completely absorbed within the top-layer. The dust particles 
re-radiate the energy at IR wavelengths, so half the redistributed 
energy leaves the disc into space, and the other half penetrates the 
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mid-layer and heats it. As a third distinct region we define the disc 
photosphere which encloses the uppermost part of the top-layer and 
whose location is defined by the parameter h (see below). 

Further, we make the following assumptions: the disc is in 
hydrostatic equilibrium in z-direction, gas and dust are well mixed 
throughout the disc, and in the mid-layer gas (g) and dust (d) are 
in thermal equilibrium at the same temperature T g = Td = Trnid- 
The emission and absorption of the gas component is ignored. 

In addition, our (gas) disc extends down to the stellar surface 
with no inner hole. The existence of such a hole would lead to 
a puffed-up inner rim and self-shadowing effects of the disc (see 
Dullemond, Dominik & Natta 2001 and Dullemond 2002). 

For a razor-thin passive disc (see Fig. 10) around a star with 
effective temperature T* and radius 7?, the monochromatic flux 
entering the disc perpendicular to the surface at distance r from the 
star is 



B„(T t ) sin 2 8 sin <i>d0 d<t> 



= B„{T,) 



R*_R 1 \_ (R*\ 
r r V \ r J 



with 




(1) 



(2) 



We assume that the star emits a Planck spectrum. At distances large 
compared with the stellar radius, this relation reduces to 



(3) 



The flux can be parametrized in terms of the solid angle fi under 
which the star is seen from a disc surface element at distance r and 
the so-called grazing angle a gr , i.e. the mean angle of incidence of 
the stellar flux 



The grazing angle of a razor-thin disc is therefore 

' R, \ R, 



(4) 



(5) 



which for large distances reduces to the handy formula first intro- 
duced by Chiang & Goldreich (1997) 



37T r 



0.4 



(6) 



This relation also holds for a wedge-shaped disc as long as its open- 
ing angle is small. 

In a flared disc the surface is curved and the grazing angle 
increases with distance from the star (see the solid lines in Fig. 
11). Therefore, the flared disc intercepts more stellar light and is 
heated more efficiently. The grazing angle of the flared disc is by 
the amount a^ rc higher than for a flat disc (Fig. 11). This flaring 
term of the grazing angle can be computed from 



flare , dh 

a„ = arctan — — 
dr 



arctan ■ 



dr 



and the total grazing angle becomes simply 



razor . flare 
■ Q'tir I ^srr 



(7) 



(8) 



disk 



i Y 



Figure 10. Sketch of stellar illumination of a razor thin disc. The grazing 
angle ct gr (see text for details) is also indicated. 




Figure 11. Sketch of a flared disc. The flaring term of the grazing angle is 
indicated. The inner radius of the dust disc is ro . 



With this new grazing angle, which is a function of h(r), i.e. the 
bottom of the photosphere or half the thickness of the disc, we can 
calculate the flux penetrating the flared disc in ^-direction. 

The photosphere encloses the uppermost part of the top-layer. 
Its onset is described by the parameter h, which is defined as the 
height z above the mid-plane where the visual optical depth along 
the direction of the infalling stellar light, i.e. along the grazing an- 
gle, equals 1. This leads to a vertical optical depth of 



Tv 



(9) 



Since the disc is assumed to be in hydrostatic equilibrium, r v can 
be calculated 



tv = KVPo(r) J e 2H^ dz (10) 

h 

where po is the density in the mid-plane and H is the scale height 



of the disc given by 



H = 



kT 



3/2 



(ID 



From equaling Eqs. {9) and l |10t h can be determined, but it is a 
function of temperature and grazing angle. 

The downwards and upwards directed fluxes, and F\ are 
in equilibrium throughout the passive disc. We can calculate these 
two fluxes explicitely at the boundary between the isothermal mid- 
layer and the top-layer. All infalling stellar radiation is reprocessed 
within the top-layer. We assume that half the incident stellar flux is 
re-radiated into space and the other half penetrates the mid-layer. 
The downwards directed flux becomes 



and the flux leaving the mid-layer in upward direction is 



(12) 
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/ j B u {T mid ) (1 



(13) 



where T m ;d is the isothermal temperature of the mid-layer, and we 
set (i = cos 9 with 9 as the angle measured from the z-axis. The 
visual optical depth of the mid-layer in vertical direction can be 
written in the form 



mid / \ mid / \ I 

t v (r) = t v (r ) — 



(14) 



with the visual optical depth at the inner edge, T™ (ro), and the 
exponent s as free parameters, and the optical depth at frequency v 
is simply 



(15) 



where <r cxt is the dust extinction cross section. The temperature of 
the mid-layer, T m id, follows from equaling Eqs. d 12b and d 1 3 b - This 
temperature is, however, a function of the grazing angle. 

We have now three important parameters, a gr , h, and T m id, 
but none of them can be computed independently, instead they must 
be calculated iteratively. 

We have to specify the visual optical depth in the top-layer in 
z-direction, T v op . It should be high enough so that first, the infalling 
stellar radiation is completely absorbed and second, the dust grains 
at the bottom of the top-layer reach a temperature close to that 
of the mid-layer in order to guarantee a smooth transition. On the 
other hand, r v op must be small enough to allow the dust emission, 
which occurs at infrared wavelengths, to escape the top-layer. 

The radiation transfer equation of a plane-parallel slab is 



I v (ji, Tv) =Ioe 



S v {t)e 



-(t— T„)//i 



dt 



(16) 



with the source function S v , the incident intensity Jo and fi = 
cos 9. The radiation transfer equation is split into up-streams, J„ , 
that penetrate the top-layer from the mid-layer, and down-streams, 
I~ ', that cross the top-layer starting from the surface. The inci- 
dent intensity is either the stellar radiation for downwards directed 
streams (but only for angles 9 under which the star can be seen, 
else it is zero), or the emission of the mid-layer, -B^(T m id)(l — 
e _T " for the upwards directed streams. 



With the help of the Feautrier parameters 



the mean intensity J„ becomes 
J v = J v {r v ) = I u v (fi,T„)dn 



(17) 
(18) 

(19) 



which is needed to calculate the source function and the emission of 
the grains. The source function itself determines the up- and down- 
streams of the intensity which leads to an iterative calculation of 
the radiation field. 

Finally, the temperature of the grains at each location in the 
top-layer follows from balancing their absorption from the sur- 
rounding radiation field, J„, and their emission at their equilibrium 
temperature, Td. 

In our code, we are adopting a mixture of spherical silicates 
and amorphous carbon grains which follow the MRN grain size 
distribution (Mathis, Rumpl & Nordsieck 1977). The particle sizes 



range from a m m = 75 A to o m ax — 1 A™. The absorption ef- 
ficiencies of the dust particles are calculated with the Mie theory 
for spherical particles and scattering of the stellar light within the 
top-layer is included. 

For the modelling of CD-42° 11721, assuming a flared disc 
scenario, we have used the physical parameters obtained by our op- 
tical analysis (extinction, effective temperature, distance, and lumi- 
nosity) and the ff-fb emission (see Sect. 6.1). We keep the top-layer 
optical depth fixed (t v ° p = 2) and use a stellar mass of 8 — IOMq 
(see Sect. 7). 

Considering this set of input parameters, the code calculates 
the location of the inner edge of the dust disc, which is defined by 
the dust evaporation temperature and may vary for different input 
values, especially of the optical depth of the mid-layer (i.e. the sur- 
face density distribution). Having fixed the inner edge of the disc, 
the disc structure and temperature distribution in vertical direction 
is calculated. Of special importance is thereby that the grain tem- 
peratures of the different grain sizes and dust components adopt to 
the isothermal temperature of the mid-layer. The dust SED is then 
calculated by integrating the emerging dust emission over the disc 
size. 

The most serious boundary condition in our calculations is the 
request of an optically thick mid-layer, justifying the application 
of a flared disc. Intending to fulfill this condition and to reproduce 
the observed dust emission, we need to play mainly with the outer 
edge of the dust disc (r out ), the mid-layer optical depth at the inner 
edge of the dust disc (r™ ld (ro)), and the exponent (s). Assuming 
values of s higher than 1.5, we obtain that only very small discs 
(f out ^ 25 AU) can keep the mid-layer optically thick, however in 
these cases, the dust disc emission is lower than the observed one. 
On the other hand, assuming s — 1 or 1.5, the density distribution 
is less steep and we can satisfy the mid-layer optical depth condi- 
tion for larger discs and consequently the dust emission is better 
reproduced. We can also get good results assuming r™ ld (?"o) in the 
range of 2.5 x 10 2 and 4 x 10 3 . Concerning the r ou t, it is limited 
by the mid-layer optically depth condition. 

The total SED is formed by the stellar, ff-fb and dust disc 
emission contributions. The stellar contribution is determined by 
the stellar parameters, especially the effective temperature and also 
the foreground extinction, obtained by our optical analysis. In the 
panel (a) of Fig. 12 we present our best result. It clearly shows de- 
viations compared to the observations, especially the near-IR part 
that cannot be well fitted by our calculations. In this model, we are 
considering a fixed density profile with s — 1, ro = 0.5AU, r out = 
60AU, and r^ id (r ) = 2.5 x 10 2 . On the other hand, an increase 
of the mid-layer optical depth (i.e. an increase in surface density) 
has caused an increase in far-IR flux and consequently, it cannot be 
adjusted to the observations. 

It is important to cite that the inclination angle of the disc, 
seen by the observer, is also an input parameter of our code and we 
have found that there is no difference in our model results, as long 
as we consider discs with inclination angles smaller than 70° (0° 
means pole-on). For larger inclination angles it is impossible to get 
a reasonable fit at all. 



6.3.2 Outflowing disc wind 

We have alternatively supposed the disc around CD-42° 11721 to 
be an outflowing disc wind, typically ascribed to B[e] supergiants, 
where the dust would be uniformly distributed in it. Our code, for 
simplification, assumes that the disc is isothermal in z-direction, 
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Figure 12. Modelling the SED of CD-42° 11721, considering in the panel 
(a) a flared disc scenario and assuming a less steep density profile (s = 1). 
As can be seen, the observations cannot be well reproduced, because the 
near-IR emission is clearly underestimated. The panel (b) show our best 
model considering the presence of an outflowing disc wind. The models 
were obtained assuming a disc seen with an inclination angle from 0° (pole- 
on) to 70° . The circles represent the photometric data obtained in the litera- 
ture. The dotted line represents the stellar contribution; the dashed line, the 
ff-fb emission; the dot-dashed line, the dust disc contribution and the solid 
line, the total emission. 



with a constant outflow velocity and a surface density described 
by: 

h h 

E(r) = J p(r,z)dz = 2 J p(r,z)dz (20) 

-h 

where h is here the height of the disc at distance r, considering a 
constant opening angle for it, and p(r, z) is the density distribution, 
being derived according to the mass continuity equation: 

M 



P(r '^ Mr 2 + «>oo 
Inserting this equation into Eq. (20), it follows that 



with 



Ma 



(21) 



(22) 



(23) 



and a is here the opening angle of the outflowing disc and ro is 
the dust sublimation radius. Eq. (22) represents the surface density 
of the gas component. To convert it into the surface density of the 
dust component, we consider the gas-to-dust ratio, that is a factor 
100. For the temperature distribution we assume a power law of the 
form 



T(r) = T 



(0 



(24) 



where To is the dust sublimation temperature, fixed at 1500 K. 

We use the same grain size distribution and grain composition 
as for the flared disc. The model output is then the integrated flux of 
the star plus disc system, with the stellar parameters taken from our 
optical analysis, and reddened with the foreground extinction. The 
inclination angle of the star-disc system is a free input parameter. 

The exponent of the temperature distribution (q, see Eq. 24) 
is found to be the parameter with most severe influence on the re- 
sulting SED. We have obtained good results assuming different val- 
ues of q in the range from to —2. For values higher than —0.4, 
we could get good results, assuming quite low disc mass loss rates 
(M disc < 10~ 7 M^yr -1 ), however the main problem in this case 
is that the temperature continues still high far away from the central 
star, what is not expected. On the other hand, using q ^ —0.5, it is 
necessary very high disc mass loss rates (M disc > lO _4 M 2/r _1 ) 
or even large disc opening angles (a > 40°) to reproduce the ob- 
served SED. 

The panel (b) of Fig. 12 shows our best result, considering 
the following parameters: a = 20°, r out ^ 20 AU, Mdisc = 
8.7 x 10~ 5 M^yr^ 1 and q — —0.4. This result seems for us better 
than those ones obtained assuming a flared disc scenario. Regard- 
ing to the outer egde of the disc (r out ), we can get the same good 
results, assuming r ou t with values from 2QAU to 8000^4(7 (which 
corresponds to approximately the size of the aperture considered). 
This means the disc is not optically thick anymore just after few 
AUs of distance from the star and the contribution from the op- 
tically thin region is negligible for the SED. This conclusion is in 
agreement with the fact that we can model the SED assuming a 
disc seen pole-on or with inclination angles lower than 70° , mean- 
ing that whole emission region is seen under these angles. Another 
interesting point is that our models are obtained assuming values 
of disc mass loss rate that are in perfect agreement with those ones 
proposed for discs of B[e] supergiants. Another parameter that also 
influences our models is the terminal velocity of the disc wind (a 
decrease of Voo causes an increase of the dust disc emission), how- 
ever based on our optical observations (see Table 1), we decide to 
keep it fixed, = v exp = 60fems _1 . 



6.4 Discussion: model results 

In the case of CD-42° 11721 several works have tentatively de- 
scribed its SED using different codes and approaches. Hillenbrand 
et al. (1992) suggested that the SED of this star shows no evidence 
of a spherically symmetric envelope. According to these authors, 
the SED could be described by a flat, optically thick massive ac- 
cretion disc. These discs present an optically thin region within 
the inner few stellar radii. These optically thin "holes" would im- 
ply the absence of warm dust which would contribute with strong 
excess of emission in the near infrared region and that is seen in 
the SED of CD-42° 11721. However they did not treat the radia- 
tive transfer problem and their fit are not good, even at 3 pm. Natta 
et al. (1993) calculated models to CD-42° 11721 taking the radia- 
tive transfer into account. Their model was solved considering a 
spherically symmetric dust envelope which would reprocess the ra- 
diation of the central source. The difference in this work is that 
they adopted a central source composed by a star and a circumstel- 
lar disc. In addition, they modified the dust properties at A > 100 
pm in order to match the observed submillimeter fluxes. Anyway, 
the model fits relatively well the data, except at longer wavelengths 
than 100 pm. 

Radiative transfer models were also calculated by Henning et 
al. (1994) with a spherical geometry for the dust envelope. Their 
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model fits well the overall SED, including longer wavelengths 
than 100 ^im. They suggested that the main part of the radiation 
comes from a spherical dust configuration. More recently, Elia et 
al. (2004) modeled the continuum emission from the circumstellar 
environment of some sources using the same code as Henning et 
al. (1994). Their models were fitted to ISO data, considering both, 
gas and dust arranged in circumstellar spherical shells with temper- 
ature and density profiles described by power laws, T ~ r~ q and n 
~ r~ p , respectively. They also considered a disc geommetry but 
they did not fit longer wavelengths. 

In our study we have firstly analyzed the photometric data 
available in the literature and based on the aperture sizes used, 
we have selected only those points that could be related to star 
(or the close-by regions) and not contaminated by other sources 
and/or nebula. Using our code which treats the radiative transfer via 
the Monte Carlo method for a spherical circumstellar medium, we 
could not to reproduce the SED of CD-42° 1 1721. Even using some 
physical parameters cited by Henning et al. (1994), especially the 
envelope inner radius, the effective temperature, and density law (n 
~r~ ) we have obtained unrealistic high grain temperatures (about 
7000 K), meaning that the dust could not be located so close to the 
star. We are aware of the fact that there are several differences be- 
tween both codes, however we can point out three of them. Firstly, 
their model is dependent of the distance of the object (luminosity) 
and they are assuming a higher value than us for CD-42° 11721 
(d = 2kpc), and secondly they considered two different density 
distributions inside of the envelope, working at different distances 
from the star. However not only the structure of the circumstellar 
region determines the shape of a SED, but also the nature of the 
grains itself. Henning et al. (1994) adopted a fluffy dust model to 
describe dust particles in such a circumstellar medium. No small 
particles were considered. That is the third difference between our 
model and Henning's one, and it may be the most important. Usu- 
ally a fluffy dust model is used to describe the coagulation process 
in dense cores of molecular clouds (Ossenkopf 1993) and young 
pre-main-sequence stars, but it is not applicable to hot stars. There- 
fore, we conclude that a spherically symmetric scenario is ruled 
out. 

Based on our optical analysis and the results described previ- 
ously, we believe that a non-spherical circumstellar envelope, pos- 
sibly an equatorial disc, seems to be necessary to explain the SED 
of the CD-42° 11721. Considering a flared disc scenario, we can- 
not model well its SED, especially the near- and the far-IR regions. 
On the other hand, considering an outflowing disc wind, we could 
get better results. It is important to cite that a better and even more 
conclusive modelling will be obtained only when more data using 
small apertures, especially from far-IR region, are available. 



7 THE NATURE OF CD-42° 11721 

As we presented in this paper, CD-42° 1 1721 is a very curious ob- 
ject presenting the B[e] phenomenon. Its evolutionary stage is quite 
doubtful, being considered either a young pre-main sequence star 
(HAeB[e]) or an evolved post-main sequence supergiant (B[e] su- 
pergiant) and all this doubt is caused by the completely absence of 
reliable physical parameters for this star. 

In our work, we have firstly presented an optical analysis of 
this object, that is based on new low- and high-resolution (highest 
than previously published) optical spectra of CD-42° 11721. Espe- 
cially the high-resolution spectrum reveals a huge number of emis- 
sion lines from both, permitted and forbidden transitions. Many of 
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Figure 13. Top panel: The HR diagram for B[e] supergiants in the Large 
Magellanic Cloud (filled circles), the Small Magellanic Cloud (open cir- 
cles), and the Milky Way (asterisks), data from Lamers et al. (1998) and 
referen ces therein. Also plotted are evolutionary tracks from Scha ller et alj 
119921) for non-rotating stars with solar metallicity. The location of CD- 
42° 11721 is in agreement of being a post-main sequence object having 
evolved from a progenitor star with an initial mass of about 8-10 Mq. 
In addition, it lies in the region of well known B[e] supergiants. Bottom 
panel : The HR diagram for galactic HAeBe stars (filled circles, Berrilli et al] 
1 19921) - The solid line defines the zero-age main sequence, the thick solid line 
defines the birthline. CD-42° 11721 is plotted as the triangle. It is located 
above the birthline. 



these lines have been detected for the first time. We created a spec- 
tral atlas containing the equivalent widths of all measured lines, as 
well as a profile classification. Line profiles are found to be single- 
peaked, double-peaked, or (especially lines of Fell) even multiple- 
peaked. The presence of atoms mostly neutral or singly ionized 
and the variety in line profiles speak in favour of a non-spherical 
density distribution within the circumstellar medium close to the 
central star. The existence of a circumstellar disc is suggested es- 
pecially by the existence of numerous emission lines from neutral 
elements. 

Using these optical data, we could estimate the stellar pa- 
rameters of CD-42° 11721. We plotted the effective temperature 
and luminosity obtained by us, into the HR diagram (top panel 
of Fig.|13l> together with known galactic and MC B[e] supergiants, 
taken from the summary table given in Lamers et al. (1998) and ref- 
erences therein. Also included are evolutionary tracks taken from 
ISchaller et all (l992) for non-rotating stars with solar metallicity. 
The position of CD-42° 11721 falls into the region of B[e] super- 
giants (L* ^ 10 4 Lq, Lamers et al. 1998) and even coincides with 
the position of Hen S59, a LMC B[e] supergiant. While all to date 
known galactic B[e] supergiants have i, ) 3 x 1O 4 L0, there 
are some low luminosity objects (10 4 < L*(L Q ) sC 3 x 10 4 ) 
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in the Magellanic Clouds (Gummersbach, Zickgraf & Wolf 1995), 
for which CD-42° 1 1721 might be the first galactic low luminosity 
counterpart, having evolved from a progenitor star with an initial 
mass of 8-10 Mq. With such a mass, the surface gravity becomes 
log g ~ 3. 

To test a possible pre-main sequence nature, we plotted the 
HR d iagram for galactic HAeBe stars taken from iBerrilli et al.l 
dl992h together with the birthline, the zero-a ge main sequence, and 
pre-m ain sequence evolutionary tracks from Bernasconi & MaedeJ 
dl996h . While the galactic HAeBe stars lie well below the birthline, 
this is not the case for CD-42° 1 1721, making its classification as a 
HAeBe star questionable (bottom panel of Fig.|T3t. 

Intending to improve this discussion, we have also made an 
IR analysis of this star. As shown in the Sect. 4, Acke & van den 
Ancker (2006) showed that the most detached characteristic of the 
ISO and IRAS spectra, the very intense PAH features, actually do 
not seem to be related to the star itself, but to an envelope or a 
surrounding nebula. This fact is very important, because the pres- 
ence of C-rich material implies in a C-rich photosphere (Waters et 
al., 1998) that seems hard to be explained in a supergiant scenario 
whose spectrum is dominated by several H recombination lines, 
even knowing that PAH emission was already identified in LBV 
stars (Voors 1999). In the Sect. 6, we have also presented that our 
best results for the modelling of CD-42° 1 1721 were obtained con- 
sidering an outflowing disc wind, typical of B[e] supergiant stars. 

Summarizing, based on our optical and IR analysis (even con- 
sidering that our errors could be underestimated) and the modelling 
of the SED, a B[e] supergiant nature seems to us to be a good 
scenario to explain the observed characteristics of CD-42° 11721. 
However, a young or even pre-main sequence nature for this curious 
object could not definitely be excluded by our qualitative analysis. 
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APPENDIX A: THE OPTICAL LINE LIST OF 
CD-42° 11721 
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Table Al. Catalogue of lines identified in the FEROS spectrum of CD- 
42 11721. Listed are the observed wavelength, equivalent width, descrip- 
tion of the line profile (S = single-peaked, D = double-peaked, M = multiple- 
peaked, A = absorption, blend = superposition of several lines), and the line 
identification with the laboratory wavelength. 



Wavelength (A) 


W(A) 


Line Profile 


Identification 


3835.3 


1.64 


s 


H9 3835 4 

117 JOJ J.^t 


JO ~J J -<J 


1 46 




Jl 11 ^1111 J JOJU.U 


3889.0 


3.56 




He I tm2} 3R88 7 

J. l\/ 1 1 111- J JOOO. / 








HR ^SSQ 1 
no jou7. i 


3933.5 


0.34 


A 


Ca II tm1 ^ 3933 7 

V. lI 11 ^111 1 / J7J J. ' 


3968.4 


13 




V lI 11 ^111 1 J J7UO.J 


3969.4 


2 33 




He 3970 1 


4068 3 


22 




[O 11J ^lllll^^ HVJUO.U 


4101.7 


4.61 




HA 41 01 7 


4173.3 


0.51 


s 


Fe II fm27"l 4173 5 

1 Gil 1 1 1 1 _ / J *"T 1 1 JJm^J 


4178.5 


0.75 


s 


Fell fm211 AMI 1 

1 t II (111^1 1 *Tl / / ./ 


4233.0 


1.70 




Fe II fm271 4233 2 

1 yZ 11 V 1 1 J / J *T^__J J 


4243.4 


0.50 


s 


[Fe II] (m21F) 4244.0 


4258.0 


0.22 


s 


Fe II (m1R\ 4258 2 

1 Gil ^lll^OJ ""T^J (J . i. 


4276.4 


0.14 




IFe III fm21 F1 4276 8 

|_1 <j 1 1 J i ] 1 1 _ 11 J a I yf . o 


4287.0 


0.57 




TFe III fm7F1 4287 4 


4290. 1 


0.13 




Ti II fm41 ^ 4290 2 

±1 11 1 111^ 1 / ^ jLy\J . jL 








Ne II fm57^ 4790 4 

lit 11 1J11.J / J jLy\J 


4296 3 


36 




Fp TT (m?&~) 47Q6 6 

1 C 11 ^lllZ.O J ti^U.U 


4303.0 


0.48 




Fe II (mil) 4303 2 

1 yz 11 V 1 1 J / f i~J\J~j 


4340.5 


8.76 




4340 5 


4351.5 


1.24 


s 


Fe II (m7T\ 4351 8 

1 yZ 11 1 111— / J tJJ 1 . (j 


4358 8 


45 




TFpTTI fm7F1 4^Q ^ 
[i ciij ^in/i j tjj y .j 


4384 7 


131 




Mtr IT 4^84 f, 
IVlii 11 tJOi.V) 


4413.5 


0.62 




TFe III (m7F1 441 3 8 

|_1 t^llj ^111/1 J 11 1 J . O 


4416.4 


2.15 


g 


TFpTTI fm^F'l 4414 "S 

[1 C 11 J ^lllUl ) 11 li'.J 








TFp TTl (m6F~) 441 ft ^ 

[1 C 11 J ^lllUl J 't't 1 U.J 








Fe II <m2T) 441 6 8 

1 r; H yill^ / J 11 1 yf . O 


4451.7 


30 




TFp TTl ('m^?F'l 44^7 1 
Li c iij ^iiuz.i j 7tJi. i 


4457 7 


18 




TFp Til (mf\V\ 44SS 


4472.4 


0.71 


s 


HeT I'm 14^ 4471 7 

11C 1 lllll^T/ I I / 1 . / 


448 1 1 


98 




Mo TT I"m4~l 44S 1 ^ 

IVlg 11 ^lll*T ) ttu 1 .J 


4489 


50 




TFpTTI ^m^F"! 4488 % 

[1 C 11J ^lllUl } ttOO.O 


4491 2 






FpTT ^m^7~l 44Q 1 4 
JT^eil ^IIU / ) L \ L ry 1 . L \ 


4508 1 


80 




Fp TT f'm^S'l 4S0R ^ 
JT^e 11 ^IllJOj 4JUO.J 








Fp TT Tm^Fl 4SO0 & 

l C 11 [111U1 J "JV77.U 


4515.1 


1.01 


M 


TFelll tm6F1 4514 9 

|_1 C11J ^lllVJl J T^JlT.7 


4520 


1 05 


D 


Fp IT 4S90 ? 

1 C 11 111JJ / J ^tJit/.Z, 


4522 4 


1.08 


D 


Fp TT ^m^S'l 4S9? ft 


4533 9 


25 




Fp IT 4S^4 ? 

1 C 11 111JJ / ) TjJi'.i 


4541 3 


36 




Fp IT ^m'^S'l 4S41 S 

1 C 11 ^lllJO } lJ4 1 .J 


4549 3 


2 07 


D 


Fp IT ^m'^S'l 4S4Q S 

1 C 11 l^lllJO } HJt7. ~j 


4555.6 


1.14 


M 


Fe II fm37"l 4555 9 

i yz ii \^iiu / j *tjjj.7 


4558.3 


0.36 


D 


Fell(m20) 4558.6 


4576.1 


0.29 


D 


Uid 


4583.5 


2.71 


blend 


Fell(m37) 4582.8 








Fell(m38) 4583.8 








Fe II (m26) 4584.0 


4588.0 


0.18 


S 


Crll(m44) 4588.2 








Crll(ml6) 4588.4 


4618.5 


0.17 


S 


C II 4618.9 


4620.3 


0.32 


s 


Fe II (m38) 4620.5 








Nll(m5) 4621.4 


4629.1 


2.03 


D 


Fe II (m37) 4629.3 


4634.5 


0.37 


s 


Fe II (m25) 4634.6 


4639.4 


0.11 


s 


On (ml) 4638.9 


4657.0 


0.24 


s 


Fe II (m43) 4657.0 


4666.5 


0.39 


D 


Fell(m26) 4665.8 


4727.7 


0.19 


s 


[Fen] (m4F) 4728.1 


4731.2 


0.40 


D 


Fell(m43) 4731.4 
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Table Al - continued 



Wavelength (A) W(A) Line Profile Identification 



4814.2 


0.28 


S 


[Fell] (m20F) 4814.6 


4823.9 


0.25 


D 


Crll(m30) 4824.1 








Sll(m52) 4824.1 


4847.9 


0.22 


S 


Fe n (m30) 4847.6 


4861.4 


30.41 


D 


H/3 4861.3 


4889.3 


0.24 


S 


[Fe II] (m4F) 4889.6 








[Fe II] (m3F) 4889.7 


4905.0 


0.16 


s 


[Fe II] (m20F) 4905.4 


4923.7 


3.63 


M 


Fe II (m42) 4923.9 


5018.1 


4.76 


D 


Fell(m42) 5018.4 


5041.5 


0.49 


S 


Si II (m5) 5041.1 


5056.1 


0.64 


S 


Si II (m5) 5056.4 


5100.5 


0.15 


s 


Fell 5101.0 


5111.1 


0.15 


s 


[Fell] (ml9F) 5111.6 


5145.6 


0.25 


s 


C II (ml6) 5145.2 








Fell (m35) 5146.1 


5149.1 


0.08 


s 


Fell 5149.5 


5158.3 


0.55 


M 


[Fell] (ml8F) 5158.0 


5163.6 


0.09 


s 


[Fell] (m35F) 5163.9 


5168.9 


3.87 


D 


Fe II (m42) 5169.0 


5183.4 


0.09 


S 


Mgl(m2) 5183.6 


5197.4 


1.98 


D 


Fen (m49) 5197.6 


5203.4 


0.10 


S 


Uid 


5227.1 


0.16 


S 


Fe III 5227.5 


5234.5 


2.41 


D 


Fe II (m49) 5234.6 


5237.0 


1.01 


S 


Uid 


5248.5 


0.25 


S 


Fe II 5248.0 


5254.6 


0.30 


D 


Fe II (m49) 5254.9 


5256.8 


0.09 


S 


[Cm] (ml3F) 5256.0 


5261.5 


0.38 


D 


[Fell] (ml9F) 5261.6 


5264.5 


0.38 


S 


Mg II 5264.3 








Fell(m48) 5264.8 


5268.5 


0.11 


S 


[Fe III] (mlF) 5270.4 


5273.0 


0.21 


s 


Cm 5272.6 


5273.1 


0.41 


s 


[Fell] (ml8F) 5273.4 


5275.7 


1.85 


D 


Fe II (m49) 5276.0 


5283.8 


0.57 


D 


[Fell] (m35F) 5283.1 


5291.6 


0.10 


S 


Fe III 5291.8 


5296.5 


0.06 


s 


[Fell] (ml9F) 5296.8 


5298.7 


0.08 


s 


Uid 


5316.3 


3.99 


D 


Fe II (m49) 5316.6 


5325.3 


0.24 


D 


Fe II (m49) 5325.6 


5333.4 


0.20 


S 


[Fell] (ml9F) 5333.7 


5346.6 


0.21 


s 


Fe II (m49) 5346.6 








[Fell] (ml8F) 5347.7 


5362.6 


1.20 


blend 


[Fell] (ml7F) 5362.1 








Fell(m48) 5362.9 


5376.0 


0.16 


S 


[Fe II] (ml9F) 5376.5 


5402.0 


0.13 


S 


Mgn (m24) 5401.1 


5413.2 


0.21 


s 


[Fein] (mlF) 5412.0 


5425.0 


0.35 


D 


Fe II (m49) 5425.3 


5432.7 


0.22 


S 


Fe II (m55) 5432.9 








[Fe II] [(ml8F) 5433.2 


5466.6 


0.25 


s 


Fen 5466.0 


5502.8 


0.13 


s 


Uid 


5510.7 


0.11 


s 


Cm (m23) 5510.7 


5527.0 


0.19 


s 


[Fell] (ml7F) 5527.3 








[Fe II] (m34F) 5527.6 


5529.5 


0.22 


s 


Uid 


5534.5 


1.00 


D 


Fe II (m55) 5534.9 


5543.8 


0.08 


s 


Uid 


5576.7 


0.84 


s 


[Ol] (m3) 5577.3 


5746.5 


0.10 


s 


[Fe II] (m34F) 5747.0 


5754.2 


0.07 


s 


[Nil] (m3F) 5754.8 



Table Al - continued 



Wavelength (A) W(A) Line Profile Identification 



5780.4 


0.44 


A 


DIB 


5791.7 


0.08 


S 


Cl(ml8) 5793.5 


5797.0 


0.16 


A 


DIB 


5835.1 


0.07 


S 


Fe II (m58) 5835.4 


5849.7 


0.05 


A 


DIB 


5876.0 


2.00 


S 


He I (ml 1)5875.6 


5889.7 


0.65 


A 


Nal (ml) 5890.0 


5895.7 


0.58 


A 


Nal (ml) 5895.9 


5957.3 


0.50 


S 


Si II (m4) 5957.6 


5978.9 


0.85 


s 


Fe III (ml 17) 5978.9 








Si II (m4) 5979.0 


5991.0 


0.27 


M 


Fell(m46) 5991.4 


6040.3 


0.13 


S 


Uid 


6045.9 


0.21 


s 


S I (mlO) 6046.0 








O II 6046.3 








O II 6046.5 


6083.7 


0.13 


s 


Fell(m46) 6084.1 


6103.2 


0.06 


s 


Fe II (m200) 6103.5 


6125.2 


0.30 


s 


Mnn(ml3) 6125.9 


6130.1 


0.32 


s 


Mnll(ml3) 6129.0 








Fe II (m46) 6129.7 


6148.2 


0.72 


M 


Fell (m74) 6147.7 


6158.0 


0.50 


s 


Ol(mlO) 6156.0 








Ol(mlO) 6156.8 








Ol(mlO) 6158.2 


6160.8 


0.14 


s 


Fe II (ml61) 6160.8 


6172.7 


0.15 


s 


Nil (m36) 6173.4 


6195.8 


0.06 


A 


DIB 


6202.9 


0.09 


A 


DIB 


6233.2 


0.25 


M 


Fell 6233.5 


6238.3 


0.49 


D 


Fe II (m74) 6238.4 


6247.8 


1.40 


M 


Fe II (m74) 6247.6 


6291.5 


0.15 


S 


Uid 


6299.8 


1.20 


D 


[O I] (mlF) 6300.3 


6317.8 


2.17 


D 


Fell 6318.0 


6338.0 


0.13 


S 


Uid 


6347.0 


1.73 


S 


Sill(m2) 6347.1 


6356.8 


0.16 


S 


Uid 


6363.3 


0.43 


D 


[O I] (mlF) 6363.8 


6371.1 


1.39 


S 


Fell (m40) 6369.5 








Sill(m2) 6371.4 


6384.1 


2.30 


M 


Fell 6383.8 


6416.6 


0.47 


D 


Fe II (m74) 6416.9 


6432.3 


0.29 


D 


Fe II (m40) 6432.7 


6442.7 


0.49 


D 


Fe II 6443.0 


6456.0 


1.94 


blend 


O I (m9) 6456.0 








Fell(m74) 6456.4 


6465.5 


0.06 


S 


Uid 


6483.4 


0.41 


S 


Nl(m21) 6482.7 








Nl(m21) 6483.8 








Nl(m21) 6484.9 


6492.2 


1.04 


blend 


Fell 6491.3 








Fell 6493.1 


6506.1 


0.35 


D 


Fe II 6506.3 


6516.2 


1.25 


M 


Fe II (m40) 6516.1 


6547.5 


0.05 


S 


[Nil] (mlF) 6548.1 


6562.9 


194.30 


D 


Ho 6562.8 


6582.9 


0.15 


S 


[Nil] (mlF) 6583.6 


6586.4 


0.38 


D 


C I (m22) 6587.8 


6598.5 


0.22 


S 


Uid 


6613.5 


0.15 


A 


Uid 


6626.5 


0.12 


S 


Uid 


6658.3 


0.09 


S 


Uid 


6666.2 


0.21 


s 


[Nill] (m2F) 6666.8 
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Table Al - continued Table Al - continued 



Wavelength (A) 


W(A) 


Line Profile 


Identification 


Wavelength (A) 


W(A) 


Line Profile 


Identification 


6678.7 


1.08 


S 


He I (m46) 6678.2 


8109.9 


0.48 


D 


[Cm] 8110.4 


6716.0 


0.08 


S 


[S II] (m2F) 6716.4 


8113.1 


0.05 


S 


Uid 


6730.3 


0.10 


S 


[S II] (m2F) 6730.8 


8114.9 


0.11 


S 


Uid 


6776.3 


0.08 


S 


Fell (m210) 6777.3 


8215.1 


0.62 


s 


Nl (m2) 8216.8 


6814.3 


0.12 


S 


Co 6812.2 


8234.7 


0.60 


s 


Mg II (m7) 8234.6 


7001.8 


0.22 


S 


1 (m21) 7001.9 


8242.8 


0.02 


s 


N I (m2) 8242.3 








1 (m21) 7002.2 


8249.9 


0.38 


s 


P40 8250.0 


7041.9 


0.14 


S 


Uid 


8267.7 


0.17 


s 


P34 8267.9 


7065.6 


1.02 


S 


He I (mlO) 7065.2 


8292.0 


0.31 


s 


P29 8292.3 


7109.1 


0.10 


S 


Uid 


8298.6 


0.83 


s 


P28 8298.8 


7154.6 


0.47 


D 


[Fell] (ml4F) 7155.1 


8323.5 


1.13 


s 


P25 8323.4 


7171.5 


0.10 


S 


[Fen] (ml4F) 7171.9 


8333.6 


1.49 


s 


P24 8334.0 


7308.1 


0.15 


S 


Fell (m73) 7308.0 


8345.2 


2.05 


s 


P23 8346.0 


7320.1 


0.16 


s 


[Oil] (m2F) 7318.6 


8358.5 


2.70 


s 


P22 8359.0 


7329.6 


0.09 


s 


[0 11] (m2F) 7329.9 


8374.3 


2.80 


s 


P21 8374.5 


7377.3 


0.30 


s 


[Ni II] (m2F) 7377.9 


8392.2 


3.34 


s 


P20 8392.4 


7387.2 


0.33 


M 


[Fe II] (ml4F) 7388.2 


8413.1 


3.62 


s 


P19 8413.3 


7409.3 


0.34 


M 


Fe II 7409.0 


8437.7 


4.58 


s 


P18 8438.0 


7423.2 


0.13 


S 


Nl(m3) 7423.6 


8445.9 


35.62 


s 


O I (m4) 8446.8 


7442.0 


0.20 


s 


N I (m3) 7442.3 


8467.4 


6.31 


s 


P17 8467.3 


7449.0 


0.14 


D 


Fe II (m73) 7449.3 


8489.8 


1.72 


s 


Fell 8490.1 


7452.0 


0.19 


S 


[Fell] (ml4F) 7452.5 


8501.7 


9.96 


s 


P16 8502.5 


7461.9 


0.57 


M 


Fe II (m73) 7462.4 


8567.2 


0.19 


s 


Nl (m8) 8567.7 


7467.7 


0.25 


D 


N I (m3) 7468.3 


8598.2 


10.77 


s 


P14 8598.4 


7479.5 


0.07 


S 


Fe II (m72) 7479.7 


8608.9 


0.14 


s 


Fe II 8609.5 


7495.4 


0.62 


D 


Uid 


8616.3 


0.26 


s 


[Fell] (ml3F) 8617.0 


7499.8 


0.40 


D 


Uid 


8628.9 


1.08 


D 


Nl(m8) 8629.2 


7506.4 


0.21 


D 


Uid 


8636.1 


0.37 


s 


Fen 8636.6 


7514.1 


1.99 


M 


Fell(m73) 7515.9 


8648.5 


0.32 


s 


Uid 


7520.3 


0.07 


S 


Fe II 7520.7 


8655.4 


0.19 


s 


Nl (m8) 8655.9 


7533.0 


0.10 


D 


Fe II (m72) 7533.4 


8664.4 


12.00 


s 


P13 8665.0 


7561.9 


0.08 


A 


Uid 


8680.0 


1.18 


D 


Nl (ml) 8680.2 


7572.1 


0.33 


M 


Fell 7571.7 








S I (m6) 8680.5 


7578.8 


0.20 


D 


Uid 


8682.9 


0.88 


s 


Nl (ml) 8683.4 


7587.4 


0.07 


S 


Uid 


8685.4 


0.62 


s 


Nl (ml) 8686.1 


7711.2 


2.32 


D 


Fell (m73) 7711.7 


8694.6 


0.27 


s 


S I (m6) 8694.7 


7731.6 


0.43 


S 


Fen 7731.7 








Fell 8695.1 


7755.2 


0.20 


D 


Fe II 7755.6 


8702.9 


0.46 


s 


Nl (ml) 8703.2 


7764.1 


0.07 


S 


[Fe II] (m30F) 7764.7 


8711.2 


0.45 


s 


Nl (ml) 8711.7 


7773.5 


6.08 


blend 


Ol (ml) 7772.0 


8718.6 


0.52 


s 


Nl (ml) 8718.8 








Ol (ml) 7774.2 


8721.9 


0.42 


s 


Fe II 8722.4 








Ol (ml) 7775.4 


8728.1 


0.25 


s 


[C I] (m3F) 8727.4 


7780.0 


0.13 


S 


Fe II 7780.4 








Nl (ml) 8728.9 


7788.9 


0.10 


D 


Fe II 7789.3 


8750.2 


11.50 


s 


P12 8750.5 


7800.9 


0.25 


D 


Fell 7801.2 


8767.8 


0.57 


s 


Uid 


7818.2 


0.16 


D 


Fell 7817.9 


8774.4 


0.09 


s 


Uid 


7835.6 


0.18 


M 


Fell 7835.9 


8783.8 


0.32 


s 


Uid 


7851.4 


0.27 


S 


Fell 7851.9 


8805.6 


0.57 


s 


Fell 8805.1 


7866.2 


1.46 


D 


Fell 7866.5 


8813.6 


0.38 


s 


Fell 8813.4 


7877.0 


1.02 


D 


Mgll 7877.1 


8819.3 


0.35 


s 


Uid 


7896.3 


1.57 


D 


Mgll (m8) 7896.4 


8827.4 


0.98 


s 


Fell 8825.5 


7917.4 


0.60 


D 


Fell 7917.8 








Fe II 8829.8 


7970.1 


0.08 


S 


Fe II 7970.4 


8833.5 


0.26 


s 


Fe II 8834.0 


7975.5 


0.60 


D 


Fell 7975.9 


8906.1 


0.17 


s 


Uid 


7981.8 


0.11 


S 


O I 7982.0 


8911.4 


0.62 


s 


Fell 8912.5 








O I 7982.4 


8915.9 


0.29 


s 


Fell 8916.3 


8030.0 


0.18 


S 


Fe ii 8030.5 


8926.2 


2.93 


D 


Fe II 8926.7 


8083.4 


0.29 


D 


Fell 8083.9 


9014.6 


11.80 


s 


P10 9014.9 


8090.4 


0.06 


D 


S 11 (m69) 8089.9 
[Til] (ml8F) 8091.8 


9060.4 


0.72 


s 


Fe II 9060.0 
Nl(ml5) 9060.6 


8092.7 


0.14 


D 


Fe ii 8092.2 








C I (m3) 9061.5 


8096.4 


0.35 


M 


Uid 








C I (m3) 9062.5 


8101.2 


0.19 


D 


Fell 8101.5 


9070.7 


0.20 


s 


[S III] (mlF) 9069.4 


8106.4 


0.42 


M 


[Cm] (m20F) 8106.9 


9075.8 


2.11 


s 


Fe II 9075.5 
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Table Al - continued 



Wavelength (A) 


W(A) 


Line Profile 


Identification 


9094.5 


1.43 


S 


C I (m3) 9094.9 


9111.5 


0.68 


S 


Cl(m3) 9111.9 


9122.7 


2.47 


D 


Fell 9122.9 


9132.2 


1.40 


S 


Fell 9132.4 


9177.1 


3.20 


D 


Fell 9178.0 


9187.1 


0.98 


S 


Fell 9187.2 


9196.9 


0.76 


S 


Fell 9196.9 


9203.3 


2.89 


s 


Fell 9203.1 
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